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Oh 1 Thermal field theory is reviewed briefly. It is noted that until recently it was not 

', known what type of real-time Green function was being calculated in the Euclidean 

ly-^ ■ approach. The formal answer to this question is then given and the unexpectedly 

T — I \ complicated answer is discussed. The physical implications of the results are then 

considered. 
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CNJ ■ 1. Introduction 

-xt- : 

Thermal field theory, or finite temperature field theory as it is also called, has been 
Q\ ■ around for at least forty years9. The resulting tools, Feynman rules etc., are very similar to 
, those found in 'normal' (zero temperature) quantum field theory, and which has been used 
Ph! very successfully in particle physics and elsewhere. However thermal field theory encodes 
o additional physics not present in zero temperature field theory, and the two are quite different 
in many ways. 

Thus, the key message of this work is that one should not be seduced into thinking thermal 
field theory is 'just' like ordinary field theory. Hence one must think carefully about the link 
between a calculation of Green functions etc. and a given physical problem and one must not 
rely upon the superficial similarities between the zero-temperature and thermal cases. 



1.1. What is thermal field theory? 

Thermal field theory is the combination of two theories. The first element is quantum 
field theory which is used to describe the behaviour of small, often fundamental, particles. 
One uses it to describe the dynamics of just of few such particles e.g. as in e + -e~ collisions 
at LEP. A fundamental aspect of the physics encoded in this theory is the idea of quantum 
fluctuations, as characterised by Heisenberg's uncertainty relation. The effects of such quan- 
tum fluctuations on physical processes is represented by Feynman diagrams with internal 
lines such as in figure |l|. The second theory is thermodynamics/statistical mechanics, which 
is used when one wants to study many body problems. This involves the ideas of statistical 
fluctuations, as seen in the appearance of Bose-Einstein, rib, or Fermi-Dirac, rif, distributions. 

The combination of quantum field theory and statistical mechanics is called thermal field 
theory. It is therefore used in many-body problems where the underlying dynamics is de- 
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scribed by quantum fields. It is describing quite a different physical situation as compared to 
that appropriate for normal quantum field theory. Both quantum and statistical fluctuations 
must be accounted for at the same time. It is therefore not surprising if thermal field theory 
has some quite different properties from standard quantum field theory. What is amazing is 
that formally the two can be described in a very similar way, as I will briefly outline below. 

As an example of the different physics, consider an electron propagating through a vacuum 
and through a QED plasma. The usual quantum fluctuations are always present and must be 
accounted for. The first Feynman diagram of figure |l]is one such contribution and the internal 
lines represent interactions of the electron with virtual particles present in the vacuum. For 
the many-body problem only, one must also consider interactions with the real physical 
particles present in the plasma. The last diagram of figure [I] represents such an example, 
while the middle two diagrams represent the effects of a mixture of quantum and statistical 
fluctuations. The graphical representation of figure |l] is closely linked to terms seen when 
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Fig. 1. Different types of process which occur when in a heat bath. 



the Feynman integral representation of the electronself-energy for Minkowskii time-ordered 
Green functions in thermal field theory is examined^. 

1.2. An example of the problem 

Consider the example of a decay rate of a particle in a many-body situation, say a baryon 
number violating process in the early universe. Typically in a one-loop Feynman diagram 
calculation, the decay rate at high temperatures is found to be of order g 2 T where g is the 
gauge coupling and T the temperature. 

However we know from our experience with zero-temperature field theory that the coupling 
we measure at lab scales is not always a good measure of the coupling strength relevant 
for processes at other scales. One excellent solution is to use renormalisation group which 
gives a link between the best values to use for the coupling constant at different scales. 
It is therefore vital to supplement the original perturbative calculation with what is non- 
perturbative information obtained from renormalisation group analysis. 

In a many-body problem, it is therefore essential to use the renormalisation group to 'run' 
the coupling from the known lab (zero temperature) value, to a value appropriate for a hot 
plasma, effectively g — ► g(T). This can be done in principle and one needs to know how 
the quantum and statistical corrections to the coupling. This usually involves calculating 
one-loop three-point diagrams. 

Now the problems begin to appear. The calculations have been done for QCD but the 
results are generic. There are two distinct calculational schemes used in thermal field theory, 
RTF (the Real-Time Formalism) and ITF (the Imaginary Time Formalism). In this case 
the calculations using RTF suggested that such diagrams could grow as fast as T 3 at high 
temperaturesa for the case where the external legs of the diagram were at zero energy. On 
the other hand, it was known that if one used ITF to do the calculation, with zero energy on 
external legs, the diagrams^! can not grow faster than T. In solving this apparent disparity, a 
proper understanding of the ITF calculations in particular emerges. 



2. A general approach to thermal field theory 



In order to provide a precise answer to the question of what is being calculated in the 
different thermal field theory schemes, I will first outline what these schemes are. I will use 
a path ordered approach to thermal field theoryUii as opposed to the alternative operator 
style used in Thermo Field DynamicsBi. The results hold in either case. The great advantage 
of this path ordered approach is that we can obtain different formalisms through different 
choices of a curve C in the complex time plane. The derivation will show that they must 
contain the same physics, they merely encode it in different ways. 

The starting point is the generating functional Z[j] 

Z[j] := Tr{e- pH T c exp[i J^dr J d 3 x j(r, x)0(r, £)]} (1) 

for the example of a single scalar field. The T c indicates some sort of time ordering of the 
operators, and the time integration is along some path C, both to be discussed later. Note 
that the statistical mechanics element is suggesting the Tr{e _/3H . . .}, so that Z[j = 0] is the 
partition function of statistical mechanics. The quantum field theory comes in through the 
path ordered exponential. This means that the Green functions, the key tools of Quantum 
field theory, are generated by taking functional derivatives with respect to the sources of the 
generating functional and then setting these unphysical sources to zero 
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Now the problem in trying to link the TFT expression ([!]) to those normally encountered 
in regular QFT is the thermal weight e~^ H . The way this is formally accounted for is the 
key trick of thermal field theory, and that is to regard the thermal weight as a time evolution 
operator. We have 

Tr{e~ pH ...} = ^{n;T in \e~ m ...\n;T in ) =^{n]T in - ...\n;r in ). (3) 

n n 

Thus we are treating this thermal weight as if it was a time translation of amount (3 in the 
Euclidean direction. Note that temperature has not replaced time in any sense, the dots in 
(IH) represent combinations of field operators which still have any time lying on the path C. 

This can now be compared with the expressions usually encountered in quantum field 
theory which are of the form (n out ; r out \ . . . \n in ; r in ). One can then just use the standard tricks 
of quantum field theory to proceed. For instance, the path integral can be used and this gives 

z b\ = J2(^iri(x);T in - ip\exp{T c J^dr J d 3 x j(T,x)(f)(T,x)}\$ in (x);T in ) 

D<f> expiz ( dr I d 3 x{£ + ]<5>)} (4) 

PBC JC JV 

where C must now run from r in to r out = T in —i(3. The PBC on the path integral stands for the 
periodic boundary conditions which must be enforced, namely that r in ) = r in — i(3) = 

Equation (|]) is a remarkably compact result. It is also clear that temperature only appears 
as a boundary condition, viz we only consider fields configurations which are periodic in time 
over —i/3. This does not stop us from having field configurations with any sort of time 



dependence in the real physical time direction or over the unphysical range of imaginary 
times up to — i/3. 

For the purposes of the problem at hand, the most important point is that all times must 
lie on a curve C which runs from Tj n to r out = Ti n — i/3. The fields in this formalism, e.g. 
the operators in Green functions generated from Z\j], only take times whose values lie on C. 
Further, working from the definition of the path integral, we see that the path ordering in ([!]) 
must be path ordering with respect to this same curve C. Thus the Green functions generated 
from Z[j] @ are ones where the field operators are ordered according to the position of their 
time arguments on C, the closer to r in (r ou t) the closer that field is placed to the left (right) 
in the thermal expectation value. Thus we have identified for the general case exactly what 
is being calculated in thermal field theories. 

The great advantage of this path ordered approach to thermal field theory is that we are 
free to choose C as we wish so long as it runs from r in to r out = r in — i/3. Thus we can 
obtain the two main formalisms through different choices of C. The derivation shows that all 
formalisms must contain the same physics, they merely encode it in different ways as they 
generate different types of Green function, path ordered with respect to different curves. 

3. The Real-Time Formalism 

The RTF (real-time formalism) can be generated by choosing the curvei of figure 0, and 
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Fig. 2. Curves used in RTF. 



then taking the limit of t{ n — > oo to give the second version in figure The curve of figure || 
curve avoids various problems found0 when alternative curves 0i are used. The parameter 
a is unphysical and does not appear in any physical quantity. 

The RTF is designed to be very like standard Minkowskii zero temperature field theory. 
When a field lies on the C\ section it is just like those encountered in Minkowskii field theory. 
Even when the field lies on C2, the field still depends only on a real-time parameter which 
varies between —00 and +00. The formalism contains fields at real physical times at the 
expense of having to give every physical field (the bits associated with G\ time arguments, 
'type 1') a thermal partner (associated with the C 2 section, 'type 2'). With experience, it is 
easy to cope with this doubling. 

The important element for the problem at hand is that this means that if we look at 
diagrams with type 1 external legs this corresponds to calculating 

Ti{e^ H T c Mri)Mr2)-Mr N )} = ^{e-^T^^)^^)...^^)}, (5) 

where = G 3?e. RTF was intended to look like zero temperature field theory, so it is 
not surprising to see that it is calculating time-ordered Green functions of fields with real 
physical times. This is what is calculated in most RTF calculations. 



4. The Imaginary-Time Formalism 

It turns out that the older ITF§ffl (Imaginary-Time Formalism) is the one which had 
not been studied in depth. The path chosen for ITF in path ordered approach just runs down 
the imaginary-time axis as shown in figure |3|. This means that only Euclidean times appear 
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Fig. 3. Path used for the Imaginary Time Formalism. 

directly in this formalism. This is fine for time independent quantities such as specific heats 
and other thermodynamic things. 

However, when one wants to learn about real dynamics we need physical, that is real, 
times. ITF, however, gives N-point Green functions which are functions of iV— 1 independent 
Euclidean times. Just like any other function one may attempt to analytically continue it 
to real values. This is usually straightforward given an explicit expression. One question 
that immediately comes to mind is what type of Green function does one obtain after such 
a continuation? After all one starts with a Green function which is the expectation value 
of fields ordered according to their Euclidean time arguments. In what order are the fields 
when all their values have been continued to real times? The answer is surprising and not 
very obvious. 

4-1. The analytic continuation of Euclidean Green functions 

Such a fundamental question for ITF has been studied for the most important case, that 
of the two-point function. This was done by Baym and Mermin in 196li'i. There they showed 
that one obtains the retarded/ advanced two-point Green functions. Further it is easy to see 
that the real part of these Green functions, which is related to the effective thermal mass, 
is the same as the real part of the time-ordered Green function. The two different types of 
Green functions do differ by 0(T) but only in their imaginary part. 

What is surprising is that this work had not been extended. The example of the running 
coupling constant, given in section 1.2, shows that understanding higher Green functions is 
important. It was as if it was assumed that the same results (real parts the same, imaginary 
parts slightly different) held for higher Green functions. 

It is not trivial to generalise the work of Baym and Mermin to N-point functions. Further, 
the results for two-point thermal Green functions do not generalise in a very obvious way 
either. However their basic approach still works. One uses three ideas. First the definitions 
of the various types of Green function are written down for a general case, the Euclidean 
time-ordered Green functions in particular. Next the boundary conditions are encoded in 
terms of Green functions with fields in a fixed order, the thermal Wightman functions. This 
is a generalisation of the KMSllli (Kubo-Martin-Schwinger) condition from two- to N-point 
functions. One simply uses the definition of a thermal Green function as trace, uses the 
cyclicity and then interprets the exp{— /3H} factor as a time translation in the usual way. 
Lastly, one finds that one must make some assumption about the behaviour of the Green 
functions at large real times. 



One thenobtains a spectral representation of the Green functions that are calculated in 
ITF. This iM 



2vr 



r (iV) ({4) = (ttV -1 r dh... dk N .s(h + ... + k N ). [ Pl23 ... N {{k}; 13,...) 



— oo 



I 

X- ; r X 



O2 + z 3 + . . . + z N ) - (k 2 + h + . . . + k N ) (z 3 + . . . + z N ) - (k 3 + . . . + k N ) 



X ... X 



zn — k 



N 



+ (ALL permutations of (123 ...N)) (6) 



The dependence on the three-momentum or space variables has not been written explicitly. 
The fc's are real energy variables and the z's are complex energy variables. One of the A;'s is 
clearly redundant, and an N-th z complex variable has likewise been introduced through the 
complex equation z j — 0. This has been done to make the expressions symmetric and easy 
to write down. 

All the individual details of the theory (masses, coupling constants etc.) and the approx- 
imation used to calculate the Green function are contained in the spectral densities, p. The 
dots in the arguments represent dependence on these factors. The spectral densities are the 
difference of two thermal Wightman functions, 

Pi23...iv( W; /?,-..) = T±{e- pH Mh)Mk2)Mh) ■ ■ ■ M^n)} /Tr{e-P H } 

-(-If Tr{e-? H 4> N (kN) ■ ■ ■ Mh)Mh)Mh)}/Tr{e^ H } (7) 

for the case of pure bosonic fields. For instance in the case of two-point functions one might 
finds there is just one spectral density and in simple cases this has the form (9{k)—9(—k))5(k 2 — 
u 2 ). The is a three-momentum dependent dispersion relation, y/(k 2 + m 2 ) for relativistic 
free fields. In general the dispersion relation and the spectral function's form are more 
complicatedB'El. 

The aim is to extract real dynamics and therefore to look at the real energy limit of the 
2's. However, it is clear that the expression (|6|) has cuts on the real z axes. Further, there 
are also cuts at real values of z\ + z 2 , z\ + z 2 + -23, and indeed in all possible sums of these 
complex energies. One must therefore specify from which side of the cut one is approaching 
by leaving in a small imaginary part, 

Zj -> Ej + iej (8) 

where Ej is real, ej is real and infinitesimal, and both can be either positive or negative. Note 
that it is the sign of the e/s, and the sign of all possible sums of the e/s, which tell us which 
side of each cut we are looking at. For > 4, the sign of the single e/s does not always 
specify the signs of all possible €j sumsS 

Now we are at last in a position to answer the question exactly what sort of real-time 
Green function has ITF calculated, in particular what operator ordering has we been left with 
after analytic continuation. Every distinct way of approaching all these cuts corresponds to 
a different operator ordering and so a different type of Green function. To see exactly what 
type of Green function has emerged, one merely does a Fourier transform of ([|) from E's to 
real times t, for a given analytic continuation (|8|). However, the resulting expression can be 
further manipulated by using identities such as 9(t) + 0(—t) = 1. This makes it highly non 
trivial to see if an expression is equivalent to something familiar or not. 

The resultsEl can be summarised as follows:- 



• The analytic continuation of ITF Green functions produces the thermal Generalised 
Retarded Funnpwns. Eq.s (^|) and (^) form a definition of these functions c.f. at zero 
temperatureE^'E-1 

• For the special case where one epsilon is positive, say Sj > 0, and all the others are 
negative (which fixes all possible sums of epsilons too), one obtains the j-th retarded 
function, Rj(tj] {t others}) EJEa'ta. The advanced functions are obtained when the signs 
of the epsilons are the opposite way round. These functions are discussed below. 

• No one analytic continuation of an ITF N-point Green function produces the time- 
ordered Green functionll^'lia. 

The retarded and advanced functions are very useful. They appear in the linear response 
approximation which trys to describe small perturbations to a plasmai@. The retarded func- 
tions are usually written as expectation values of multiple commutators, with theta functions. 
One field is picked out as being special, they are symmetric only in the remaining fields. That 
one special field must have the latest time if the retarded function is to be non-zero. This 
means there are N N-point retarded functions. The advanced functions are the same but 
with all the opposite theta functions 6{t) — > 9(—t) and a possible overall minus sign. For 
instance, for pure bosonic fields, 

Riih; t 2 , t 3 , U) = 0(t! - t 2 )6(t 2 - t 3 )6(t 3 - t 4 )Tr{e- /3H [[[0 1 , 2 ], 3 ], M}/Tr{e^ H } 

+ (All (234) permutations) (9) 

More recently, Baier and Niegawa have derived the link between ITF and the retarded/advanced 
functions only, in a much simpler wayta. 

One interesting sideline is that I do not know how many of these generalised functions 
ITF produces. There are N retarded and N advanced functions, but they are a small subset 
of the total number of generalised retarded functions. The results to date are shown in table 
([]]) • The results for N=2-6 were produced elsewhere^ whereas the N = 7, 8 results were 
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11,292 
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Table 1. Numbers of iV-point functions. The rows are:- N, number of external legs; R, A numbers 
of retarded and advanced functions; RTF, numbers of distinct RTF functions; TWF, the number of 
thermal Wightman functions; ITF=GRF, the number of generalised retarded functions. 

found first by van Eijck using a different method. I have cross checked the N = 7 result. The 
results up to N = 6 were independently derived in the zero temperature case by Araki0 but 
the N = Q result quoted there seems to be incorrect. 

5. Physical Implications 

The most important physical implication comes when we try to compare time-ordered 
functions and retarded/advanced functions. This is because the former are most easily pro- 
duced in RTF whereas ITF is almost always calculating the latter. One finds for trun- 
cated/lPI functions0 

r^(Si, E 2 , E 3 ) = n 2 n 3 (aii^ p/ - e~ pEl A\ PI ) + (All (123) permutations) (10) 



where rij = (exp{j3Ej} — crj)^ 1 and Oj is +1 (—1) if the j-th field is bosonic (fermionic). Aj is 
the j-th advanced thermal Green function. For the case of purely bosonic fields, we find that 
the real parts of the Green functions in the limit of T ^ Ej satisfy 

Re{T] p I (E 1 ,E 2 ,E 3 )} 

= T 2 (Re{R\ P1 }—^— + (All other (123) permutations)) + 0(1) (11) 
V E 1 E 2 E 3 J 

Roughly speaking then, we find that the time-ordered (T\ p ) and retarded/ advanced 1PI 
functions differ by a factor of the square of the temperature. This accounts for the differences 
in the results of the calculations of QCD three-point diagrams!'! described in section 1.2. 

From this work we see that the same Feynman diagram represents different types of 
thermal Green function for different formalisms. Further, different types of thermal Green 
function differ by large amounts. Therefore it is essential that one understands quite carefully 
which type of Green function is required for a given physical problem. Only then, using 
the results described here, can one choose the calculational scheme best suited for your 
purpose. The similarity in the Feynman rules of thermal and normal zero temperature field 
theory obscures the fundamental differences between the physics of the two theories. These 
differences arise because thermal field theory is modelling a very different physical situation, 
a quantum plasma not a quantum vacuum. 
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